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Abstract Quantum electrodynamic (QED) effects that shift the binding 
energies of hydrogenic energy levels have been expressed in terms of a 
semi-analytic expansion in powers of Za and ln[(Za) -2 ], where Z is 
the nuclear charge number and a is the fine-structure constant. For 
many QED effects, numerical data are available in the domain of high 
Z where the Za expansion fails. In this Letter, we demonstrate that 
it is possible, within certain limits of accuracy, to extrapolate the Za- 
expansion from the low-Z to the high-Z domain. We also review two- 
loop self-energy effects and provide an estimate for the problematic 
nonlogarithmic coefficient Bqq . 
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1 Introduction 

The subject of the current Letter is the investigation of QED radiative corrections in bound 
hydrogenlike systems which provide one of the most stringent and accurate available tests of 
quantum field theory and are amenable to high-precision spectroscopy on which the determi- 
nation of fundamental constants is based [1]. The purpose of this investigation is twofold: first, 
to demonstrate that the Za-expansion which is inherently valid only at small Z, can be ex- 
trapolated by "deferred" Pade approximants to the domain of high Z, albeit with a certain loss 
of accuracy in the theoretical predictions. The second purpose is to provide a brief review of 
logarithmic two-loop higher-order binding corrections to the Lamb shift of hydrogenic states. 
The calculation of these effects has recently been completed [2,3], but results have been provided 
only for the total effect which is the sum of the two-loop self-energy, two-loop vacuum polariza- 
tion and combined effects. For a comparison to numerical calculations which are currently being 
pursued [ ], it is helpful to analyze the coefficients that relate to specific sets of gauge-invariant 
diagrams. This second purpose is actually a prerequisite for carrying out the extrapolation of the 
two- loop self-energy from low to high Z, wherefore we commence this Letter with the endeavour 
of providing some "mini-review" of the two-loop hydrogenic energy shifts. 

2 Brief Review of Two— Loop Hydrogenic Energy Shifts 

We use natural Gaussian units with h = c = cq = 1 and e 2 = 47ra, as it is customary for 
QED bound-state calculations. The two-loop radiative shift of a hydrogenic S state, within the 
Za-expansion, reads 

= O 2 W m 

where 

H(Za) = B m + (Za) 5 5 o + (Zaf {B 63 \n 3 (Za)- 2 

+5 62 \n\Zay 2 + B 61 \n(Zay 2 + B 60 } + . . . , (2) 
and the ellipsis denotes higher-order terms. 

This section is a brief review of the known two-loop coefficients £>4o, £>5o, Bq%, Bq2 and B§\. 
The Feynman diagrams which contribute to the two-loop bound-state energy shifts are shown 
in Fig. 1. 

For S states, the contributions to B^ can be evaluated by considering the form-factor approach 
described e.g. in Sec. VIII.B.1-VIII.B.3 of [ ] or in Sec. 1 of [6]. Detailed information about the 
F{(0) form factor slope and the magnetic form factor 1*2(0) attributable to the different sets of 
diagrams (2LSE) and (SVPE) shown in Fig. 1 are given in Eqs. (16), (17) and (18) of Ref. [6]. 

We recall the following n-independent results for the £>4o-coefflcients of S states, 
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Figure 1: The Feynman diagrams that contribute to the two-loop QED 
energy shifts of hydrogenic bound states fall quite naturally into three 
separately gauge invariant categories: (i) the two-loop self-energy effects 
(2LSE), which are historically the most problematic, (ii) the vacuum- 
polarization insertion into the virtual photon line of the one-loop self- 
energy (SVPE), and (hi) diagrams involving both the self-energy and 
the one-loop vacuum polarization on the one hand and pure two-loop 
vacuum-polarization corrections on the other hand, summarized here 
as the set (SEVP). The double line denotes the bound-state electron 
propagator, i.e. including all Coulomb interactions. 



Note that the distribution of £>4o-contributions among the different sets of diagrams in Fig. 1 
[notably (2LSE) and (SVPE)] is different from the separation into a "self-energy" correction 
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[Eq. (A24) of [ ]] and a "magnetic moment contribution" [Eq. (A25) of [ ]]. We also recall that 
the first treatment of the leading two-loop self-energy coefficient B40 was completed in [7]. 

The evaluation of the relativistic correction B$q due to the two-loop self-energy has represented 
a considerable challenge [ , ]. Diagrams involving a closed fermion loop were studied in [10,11]. 
It might be useful to point out that the contribution of the diagram (SVPE) is the sum of the 
contribution labeled En and Ey in Ref. [ )]. The coefficients read: 

s£ LSE) (nS) = -24.2668(31), (4a) 

£g VPE) (nS) = -0.1571, (4b) 

B 50 EVr Vs) = 2.8677, (4c) 
B 50 (nS) = -21.5562(31). (4d) 

These results are in agreement with the data presented in Eqs. (A28) — (A29) of [ ]. 

The coefficients of sixth order in Za have recently been analyzed in [2] (see also the references 
therein). The triple logarithm Bqs originates exclusively from the set (2LSE). The diagrams of 
the set (SEVP) have been calculated in Sec. VII of [2], and the results for the double and single 
logarithms originating from these diagrams can be obtained by adding the contributions labeled 
Ey P and Ey P in Ref. [2], and the additional logarithm implicitly contained in the Eq. (40) ibid. 
which is proportional to B^^ F \ The diagram (SVPE) generates only a single logarithm given 
by B\ /2 [again consider Eq. (40) of Ref. [2]]. In total, the results read as follows: 

B 63 (nS) = #g LSE) (nS) = - A = -0.296296 . (5) 

The total result for Bq2 as well as its n-dependence were obtained in Refs. [2,12,13,14]. Here, we 
give the formulas for the particular sets of diagrams, for the case n — 1 as well as the difference 
to a state of general n: 
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where ^ denotes the logarithmic derivative of the gamma function, and C = 0.577216... is 
Euler's constant. 



4 



The formulas for Bqi are a little more involved, 
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The results for iV(nS) taken from [3] read 

iV(lS) = 17.855 672(1) 
iV(3S) = 10.449810(1) 

iV(5S) = 9.304114(1), iV(6S) = 9.031832(1), 

iV(7S) = 8.840123(1), iV(8S) = 8.697639(1) . 



AT(2S) = 12.032 209(1) 
JV(4S) = 9.722413(1), 
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The slightly shifted results for TV (IS) also explains a discrepancy in an intermediate step of the 
calculation of radiative corrections to the muonium hyperfine splitting [15,16]. 



3 Extrapolation of the Za— Expansion 

We start our consideration with the one-loop self-energy which is the dominating radiative 
correction in hydrogenlike bound systems. We write the (real part of the) one-loop self-energy 
shift A£^ L) as 

AE^ = °- {Zaf ™ F(Za) , (9) 

where F(Zol) is a dimensionless quantity which depends on the principal quantum number n, 
the total electron spin+angular momentum j and the electron orbital angular momentum /, and 
of course on the parameter Za. 
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The semi-analytic expansion of F(Za) about Za — for P states and states with higher angular 
momenta gives rise to the expression [21], 

F(Za) = A 40 + (Za) 2 [A 61 \u{Za)~ 2 + A 60 ] +... (I > 1) , (10) 

where the ellipsis again denotes omitted higher-order terms. The Aqq coefficient has proven to 
be by far the most difficult to evaluate [22,23,24,25,26], and for 2P states, results have become 
available recently [17,18]. 




Figure 2: Extrapolation of the semi- analytic Za-expansion of the one- 
loop self-energy (10) to the range of high nuclear charge via "deferred" 
[2/2]-Pade-approximants for the 2P 3 / 2 -state as described in the text. 
The analytic coefficients Aq\ and Aqq in Eq. (10) are taken from 
Refs. [17,18]. The extrapolated semi-analytic Za-expansions are closer to 
the numerical data for high Z than the "raw" Za-expansion. Numerical 
data at high Z are taken from Refs. [19,20]. 



The semi-analytic expansion (10) is generally assumed to converge to the function F(Za) for 
low Z, at least in an asymptotic sense. This is confirmed by recent numerical evaluations [27,28] 
for S and P states and the successful consistency check with available analytic results [17,26]. 
In many cases, an asymptotic expansion valid a priori for small expansion parameter Za can 
be extrapolated to large coupling, if it is combined with a suitable convergence acceleration or 
resummation method (the latter in the case of a divergent input series [29]). The logarithms in 
Eq. (10) make a power series expansion about Za — impossible. However, an extrapolation is 
still possible if we expand about Za — a ^ and use the fact that the nonperturbative function 
F{Za) 1 in the range of small is very well represented by the first terms in its asymptotic 
expansion, as suggested by Fig. 1 of Ref. [27]. The logarithms in (10), when expanded about 
Za — a, give rise to an infinite power series in the variable g — (Z — 1) a. We proceed as follows: 
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Figure 3: Extrapolation of the two-loop self-energy. Analytic data are 
taken from Sec. 2 [see the (2LSE)-parts of Eqs. (3) — (8)]. Numerical 
data are found in Ref. [4]. Much better agreement between numerical 
and analytic data is achieved for negative Bqq. 



for all radiative corrections studied in the sequel, we start from the semi-analytic expansion 
and take into account all known coefficients. We then expand in g and evaluate the diagonal 
[2/2]-Pade approximant to the resulting power series (for the definition and a comprehensive 
discussion of Pade approximants we refer to [30]). This could be characterized a "deferred" 
approximant which is evaluated only after one has "advanced" to the point Z = 1 from the 
"starting point" Z = (or equivalently Za = 0). Formally, the semi-analytic Za-expansion is 
performed about the point Z = 0. In re-expanding the perturbation series about a different 
point in the complex plane, we follow ideas outlined in [31] which were originally applied to 
the problem of calculating the autoionization width of atomic resonances in an external electric 
field. 

The choice of the [2/2] -Pade approximant is motivated by the paradigm of finding a compromise 
between the necessity to harvest the information contained in the logarithms, which give rise to 
power series terms of arbitrarily high order, and at the same time to avoid spurious singularities 
which may be incurred when the degree of the Pade approximant is increased to an excessively 
high order. The [2/ 2] -deferred Pade approximant about Z = 1 to the function F [see also Eq. (3) 
of [32]] has five parameters po, . . . ,£>2 5 qi, • • • > Q2, 

[2/2Hg)= E h Ptgl . , (11) 
1 + Ej=i Qj 9 3 

which are determined by the requirement that the power series expansion of [2/2]p(g) about 
g = reproduce the power series expansion of F(Za) with Za = a + g, also about g = 0, 
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up to the order 0(g 4 ). As discussed in [30], this condition alone defines the Pade approximant 
uniquely. The "one" in the denominator of (11) is the so-called Baker convention. 

Because we do not observe a factorialfy divergent alternating-sign pattern in the ^-expansion, 
we do not employ the delta transformation [29, Ch. 8] which has proven to be superior to Pade 
approximants in a number of applications where factorial divergence is observed (e.g., [ ]). We 
rely on the robust Pade approximants, while stressing that it may be possible to find better 
extrapolation algorithms that harvest the analytic structure of (10) and give rise to logarithmic 
terms naturally. As yet, we have been unable to find such algorithms. 

Fig. 2-3 show that the extrapolated semi-analytic expansions have a somewhat better agreement 
with medium and high-Z numerical data than the known terms of the Za-expansion alone. We 
observe that the high-Z results for the energy correction given by the irreducible set of two- 
loop self-energy insertions into the bound electron propagator (see Ref. [ ]) can only be made 
consistent with our extrapolated Za-expansion if we assume that the coefficient £>6o(lSi/2) is 
negative and rather large in magnitude. 

4 Conclusions 

We have presented a "mini-review" of recent two-loop self-energy calculations [2,3,12] in Sec. 2, 
clarifying the distribution of sixth-order (in Za) two-loop binding corrections to the Lamb shift 
over the set of diagrams shown in Fig. 1. Results for the two- loop coefficients, including excited S 
states, are provided in Eqs. (3) — (8). The distribution of the logarithmic corrections over distinct 
sets of diagrams needs to be clarified in order to allow for an accurate comparison to numerical 
calculations which are currently being pursued [4]. In Sec. 3, we present a crude extrapolation 
scheme for the extrapolation of the Zo:-expansion from low Z to high Z. The scheme follows 
ideas outlined in [ ] and is based upon an expansion in the variable g where g is defined as 
Za = a + g. The "deferred" Pade approximant is then evaluated in terms of the variable g, 
i.e. after the coupling parameter Za has acquired the value a / 0, "starting" from Za = 0. 
This deferment circumvents the problems introduced by the logarithms in Eqs. (2) and (10); 
however, we stress here that it would be highly desirable to find better extrapolation algorithms 
that harvest the analytic structure of (2) and (10) and give rise to logarithmic terms naturally. 
Although the extrapolation scheme has problems (in some cases, we observe spurious poles in 
the Pade approximant at medium-Z values), we have observed rather consistent improvement 
over the "raw" Za-expansion with this scheme for a number of states and QED effects which we 
studied using the "deferred" Pade-extrapolation scheme. Details will be presented elsewhere. 
Based on our extrapolations of the two-loop effect in Fig. 3, we would like to advance the 
tentative estimate £>6o(lSi/2) ~ —100(50). 
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